Heat Conduction in the Vortex State of NbSe 2 : 
Evidence for Multi-Band Superconductivity 



Etienne Boaknin, 1 M.A. Tanatar, 1 Johnpierre Paglione, 1 D. Hawthorn, 1 F. Ronning, 1 R.W. 
Hill, 1 M. Sutherland, 1 Louis Taillefer, 1 ' 2 Jeff Sonier, 2 ' 3 S.M. Hayden, 4 and J.W. Brill 5 

1 Department of Physics, University of Toronto, Toronto, Ontario, Canada 
2 Canadian Institute for Advanced Research, Toronto, Ontario, Canada 
3 Simon Fraser University, Department of Physics, Burnaby, Canada 
J *H. H. Wills Physics Laboratory, University of Bristol, Bristol BS8 1TL, United Kingdom 
5 Department of Physics and Astronomy, University of Kentucky, Lexington, Kentucky, 40506-0055 

(Dated: February 1, 2008) 

The thermal conductivity k of the layered s-wave superconductor NbSe2 was measured down 
to r c /100 throughout the vortex state. With increasing field, we identify two regimes: one with 
localized states at fields very near H c i and one with highly delocalized quasiparticle excitations at 
higher fields. The two associated length scales are naturally explained as multi-band superconduc- 
tivity, with distinct small and large superconducting gaps on different sheets of the Fermi surface. 
This behavior is compared to that of the multi-band superconductor Mgfb and the conventional 
superconductor VaSi. 

PACS numbers: 74.70.Ad, 74.25.Fy, 74.25. Qt, 74.25.Jb 
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Multi-band superconductivity (MBSC) is the existence 
of a superconducting gap of significantly different magni- 
tude on distinct parts (sheets) of the Fermi surface (FS). 
This unusual phenomenon has recently emerged as a pos- 
sible explanation for the anomalous properties of some 
s-wave superconductors. Although first experimentally 
observed over twenty years ago , the possible existence 
of MBSC has not often been considered since then. The 
current interest in MBSC has been fueled by the peculiar 
properties of the 40-K superconductor MgB2, where the 
case for MBSC is now rather compelling 2}. In particu- 
lar, a gap much smaller than the expected BCS gap has 
been resolved in tunneling experiments One conse- 
quence of such a small gap is the ability to easily excite 
quasiparticles, which, for example, can make the prop- 
erties of this s-wave superconductor similar to those of 
d-wave superconductors. 

Based on angle-resolved photoemission (ARPES) mea- 
surements, it has recently been proposed that the 7- 
K layered superconductor NbSe2 is also host to MBSC 
Q. A sizable difference in the magnitude of the super- 
conducting gap was found on two sets of Fermi surface 
sheets, with no detectable gap on the smallest sheet. 
However, these measurements were only performed at 
5.3 K. It is clearly of interest to shed further light on this 
sheet-dependent superconductivity by performing bulk 
measurements down to low temperatures. 

In this Letter, we report a study of heat transport in 
NbSe2 down to T c /I00 throughout the vortex state, pro- 
viding further evidence for MBSC. By measuring the de- 
gree of derealization of quasiparticle states in the vortex 
state, heat transport probes the overlap between core 
states on adjacent vortices, i.e. the size of the vortex 
core (~ £), and hence the magnitude of the gap (~ l/£)- 
We resolve two regimes of behavior: one limited to very 



low fields (up to ~ 5H c i), where derealization is slow 
and activated as in conventional (single-gap) supercon- 
ductors like VsSi, and one for all other fields up to H C 2 
where quasiparticles transport heat extremely well, as in 
unconventional superconductors with nodes in the gap. 

NbSe2, a quasi-2D metal with hexagonal symmetry, 
displays a transition to a charge density wave state 
around T ~ 35 K. Measurements of its FS at low tem- 
peratures, by both de Haas-van Alphen (dHvA) [§J and 
ARPES measurements |4| , agree with band structure cal- 
culations that predict a FS made of 4 or 5 sheets. These 
sheets divide into two groups: a small T-centered pocket 
derived from the Se 4p band (denoted as T band) and 
larger nearly two-dimensional sheets derived from Nb 4<i 
bands. Scanning tunnelling spectroscopy (STS) at 50 mK 
revealed a spectrum that is consistent with a distribution 
of gaps that range from 0.7 to 1.4 meV 6J. In the vortex 
state at very low fields (near H c x), a zero-bias conduc- 
tance peak, characteristic of states localized in the vortex 
core, was observed at the vortex center 

The thermal conductivity K of NbSe2 was measured 
in a dilution refrigerator using a standard technique @|. 
Measurements where made at temperatures increasing 
from 50 mK and in magnetic fields ranging from to 6 T, 
applied parallel to the c-axis and perpendicular to the 
in-plane heat current. The sample was cooled in field to 
ensure field homogeneity. Measurements as a function of 
field at fixed temperature resulted in nearly no difference 
as compared to the field-cooled data (see Fig. 2a). 

The sample, a rectangular parallelepiped with dimen- 
sions 1.2 x 0.5 mm in the plane, and 0.1 mm along the 
c-axis, is from the same batch as the sample used by 
Sonier et al. [lfl lllj , and has a superconducting transi- 
tion temperature T c = 7.0 K with a width ST C — 0.1 K. 
It was cleaved to provide six fresh surfaces for silver 
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FIG. 1: Thermal conductivity of NbSe2 at several applied 
fields, plotted as k/T vs T 2 . The solid line indicates the 
value expected from the Wiedemann-Franz law as obtained 
from resistivity measurements at H — 4.5 T. The field is 
applied parallel to the c-axis and perpendicular to the heat 
current Q. 



paint contacts, with resistances at low temperatures of 
roughly 20 mf2. The residual resistivity ratio is 40 
(po ~ 3 fifl cm), and the upper and lower critical fields 
are respectively H C 2 = 4.5 T and H c \ = 20 mT for H || c. 
The coherence length estimated from H c2 is £(0) = 85 A. 

The thermal conductivity of NbSe2 is plotted in Fig. 1, 
as k/T against T 2 . This enables a separation of the elec- 
tronic and the phononic thermal conductivities, since the 
asymptotic T dependence of the former as T — > is linear 
while that of the latter is cubic. The electronic thermal 
conductivity kq/T is thus obtained as the extrapolated 
T -> value. In zero field, n /T = 0.000 ± 0.005 mW K 
~ 2 cm -1 , a clear indication that NbSe 2 is an s-wave su- 
perconductor with a fully gapped excitation spectrum. 
However, by applying a small magnetic field (H > H c i), 
an electronic contribution develops as a rigid shift from 
the H = curve in Fig. 1. At higher fields (H > 1.5 T), 
the electronic contribution dominates the conduction 
over the entire temperature range and k/T is constant in 
temperature within our experimental resolution. Above 
H C 2, the Wiedemann-Franz (WF) law is satisfied and the 
thermal conductivity saturates. 
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FIG. 2: (a,b) Thermal conductivity of NbSe2 (empty circles) 
and VsSi (diamonds) at T — > vs H, normalized to values at 
H C 2- Filled circles come from a sweep in field at T — 130 mK 
from which the H = thermal conductivity (phononic con- 
tribution) has been subtracted. The thick solid line in (b)is 
a theoretical curve for the thermal conductivity of VaSi |j. 
The thin line is a guide to the eye. (c) k/T vs H for NbSe2 at 
T = 130 mK (circles) and T = 900 mK (squares). The latter 
shows a typical drop in the phononic thermal conductivity at 
H c i — 20 mT. The former shows that the electronic thermal 
conductivity starts to increase right at H c i but has a slow 
activated- like behavior for fields below 0.03 H C 2- 



The kq /T values are plotted as a function of H on a re- 
duced scale in Fig. 2. Also plotted is a field sweep at T = 
130 mK from which the zero- field value (the phononic 
contribution) has been subtracted. As seen in Fig. 2c, 
heat conduction starts to increase right at H c \ in what 
could be qualified as an activated behavior, although in 
a very limited range of fields: H c i < H < 0.03H C 2- (The 
value of H c i is determined in situ as the drop in the 
phonon k due to vortex scattering.) At higher fields, 
k /T increases rapidly, i.e. faster than (i7/i? c2 ) k^/T, 
where k^ /T is the normal state value. This shows the 
presence of highly delocalized quasiparticle states almost 
throughout the vortex state of NbSe2- 

This is in stark contrast to the behavior expected of a 
type-II s-wave superconductor. Indeed, when a field in 
excess of H c \ is applied, and vortices enter the sample, 
the conventional picture is that the induced electronic 
states are localized within the vortex cores. As one in- 
creases the field, the intervortex spacing d ~ \/&o/ B 
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FIG. 3: (a) Thermal conductivity and heat capacity of NbSe2 normalized to the normal state value vs H/H C 2- The heat capacity 
was measured in two different ways: i) at T = 2.4 K on the same crystals as used in this study fToj, and ii) extrapolated to 
T — > from various temperature sweeps on different crystals [Ts|. (b) Equivalent data for MgB2 single crystals flU |2C| . (c) 
Equivalent data for VaSi, with a theoretical curve for k/T The specific heat is measured at T = 3.5 K |2lB and extrapolated 
to T = |22^|. The straight line is a linear fit. The thermal conductivity is seen to follow the specific heat very closely for both 
NbSe2 and the multiband superconductor MgB2. It does not, however, for the conventional s-wave superconductor V3SL 



decreases. The localized states in adjacent vortices will 
have an increasing overlap leading to enhanced tunneling 
between vortices, and the formation of conduction bands. 
Strictly speaking, the electronic states are actually al- 
ways delocalized but with extremely flat bands at low 
fields 01 • As these gradually become more dispersive, 
the thermal conductivity should increase accordingly and 
grow exponentially with the ratio d/£, as is indeed ob- 
served in Nb . 

A better point of comparison for NbSe2 is VsSi, an ex- 
treme type-II superconductor with comparable supercon- 
ducting parameters (T c = 17 K, £ = 50 A). This is done 
in Fig. 2, where the thermal conductivity of VaSi is seen 
to grow much more slowly with H than that of NbSe2, as 
described by theory [9j. Quantitatively, at H = H c2 /2Q, 
pJT = i x n N /T for NbSe 2 and ^ x n N /T for V 3 Si 
[14j. Note that the samples compared in Fig. 2 are 
in the same regime of purity. From the standard rela- 
tion £(0) = 0.74£ [x(0.88£o/0] 1/2 > we obtain for V 3 Si 
io/l = 0.13, with £(0) = 50 A from H c2 and I = 1500 A 
from dHvA [J^. This is similar to the value of 0.15 for 
NbSe 2 0. 

The high level of delocalization in NbSe2 is a clear 
indication of either a gap with nodes (e.g. d-wave) or a 
nodeless gap which is either highly anisotropic or small 
on one FS and large on another. A gap with nodes is 
ruled out by the absence of a residual term in the thermal 
conductivity in zero field [l7j . 

It is revealing to compare NbSe2 to MgB2, for which 
the thermal conductivity has a similar field dependence. 
Strikingly, k follows roughly the same field dependence 
as the specific heat C for both NbSe2 and MgB2- This 
is shown in Fig. 3 where n(H) and C(H), are plotted on 
a reduced field scale for sing le crystals of NbSe2 0, 0] , 
MgB 2 [Hill, and V 3 Si HJEl, with H\\c for hexagonal 



NbSe2 and MgB2, and H\\a for cubic V3SL In conven- 
tional superconductors like VaSi, k(H) and C(H) are 
very different (see Fig. 3c) because the excited electronic 
states are largely localized. 

MgB2 is a well established case of MBSC with a small 
gap on one FS (A w = 1.8 meV) and a large gap on the 
other (A CT = 6.8 meV). The field dependence of its heat 
capacity is well understood in this context [20j, with a 
distinctive shoulder at a field of H c2 /10 (see Fig. 3b) 
|23| . A similar shoulder is also manifest in NbSe2 around 
H c2 /9 (see Fig. 3a). Empirically, the striking fact that 
heat transport and heat capacity have the same field de- 
pendence in both materials points to a common expla- 
nation, and hence sug gest s that NbSe 2 is host to multi- 
band superconductivity [24J . This is consistent with recent 
ARPES measurements at T = 0.8 T c 0. 

In conventional superconductors, the delocalization of 
vortex core bound states occurs gradually on the scale 
of H c2 , and the characteristic length scale is £(0) ~ 
V$o/27r# c2 . It appears that in NbSe2 (and MgB2), 
there are two characteristic length scales for delocaliza- 
tion: £* and £(0). To see this, we focus on the low field 
region. For NbSe2, both k/T and C/T have been mea- 
sured with high precision on the same crystals, thereby 
making a detailed comparison possible. Fig. 4a shows the 
comparison for fields below H c2 /1Q, where the two do not 
coincide: C/T increases abruptly above H c i while k/T 
grows slowly, in an activated way. This is consistent with 
the presence of localized states at very low fields as im- 
aged by STS 0, ■ Then this behavior gives way to a 
rapid increase of the thermal conductivity at fields above 
0.03 H c2 . This is a clear indication that the field scale 
associated with delocalization in NbSe2 is much smaller 
than H c2 . 

In fact, we can scale the behavior of the low- field 
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FIG. 4: (a) Thermal conductivity (circles) and the field evolu- 
tion of the heat capacity SC/T = [C(T, H) - C(T, H = 0)]/T 
(squares) [To|| of NbSe2 normalized to the normal state values 
vs H / H C 2 at very low fields. The line is a guide to the eye. (b) 
Normalized thermal conductivity vs (H — H c i)/H* for NbSe2 
and VsSi with H* = H C 2/Q and H C 2 respectively, (c) Ratio 
of heat transport to heat capacity in NbSe2. 



the r band whereas they extract an average value of 
A e -ph = 1.8 from specific heat The origin of MBSC 
may lie in the fact that A e - P h is smaller for the F band. 
This suggests that superconductivity originates in the Nb 
Ad bands and is induced onto the Se 4p T pocket. 

In summary, measurements of heat transport in the 
vortex state of NbSe2 at low temperatures reveal the ex- 
istence of highly delocalized quasiparticles down to fields 
close to H c \ . This is in striking contrast with what is ex- 
pected in a s-wave superconductor where well-separated 
vortices should support only localized states, as is ob- 
served in V3SL We identify two characteristic length 
scales that govern the destruction of superconductivity 
in NbSe2: the usual one associated with H C 2 and another 
associated with a much smaller field H* ~ H C 2/9. We 
attribute this to multi-band superconductivity, whereby 
the gap on the pocket-like T band is approximately 3 
times smaller than the gap on the other two Fermi sur- 
faces. 

We are grateful to H.-Y. Kee, Y.B. Kim, W.A. MacFar- 
lane, K. Machida, A.G. Lebed and K.B. Shamokhin for 
stimulating discussions. We thank Sasa Dukan for pro- 
viding us with theoretical calculations. This work was 
supported by the Canadian Institute for Advanced Re- 
search and funded by NSERC of Canada. 



thermal conductivity of NbSe2 to that of VaSi using 
H* = H c2 /9 (Fig. 4b). This is also seen clearly if we 
plot the ratio of the thermal conductivity to the specific 
heat (Fig. 4c) which measures the degree of delocaliza- 
tion. The ratio is seen to have two regimes: a rapid 
increase below H* ~ iJ C 2/10 and a slow one above. In 
summary, the second length scale associated with delo- 
calization in NbSe 2 is f* ~ £(0)/\/9 = £(0)/3. This is 
consistent and may explain naturally the shrinking of the 
vortex cores observed with muon spin rotation [111 ] . 

Considering the fact that the upper critical field is re- 
lated to the superconducting gap by H C 2 oc A 2 /vp where 
vp is the Fermi velocity, we estimate the gap to vary 
over the FS by a factor of 3 (A* ~ A /3). (Note that 
we assume vp to be constant, within the direction of the 
ab- plane, as found by band structure calculations 
A variation of this order was reported earlier from STS 
measurements |fj], where the spread in A was measured 
to be between 0.7 meV and 1.4 meV. Moreover, theoreti- 
cal modelling of STS features associated with the vortex 
cores j2^| and efforts to model C(T) [2(| were found to 
require an anisotropy in the gap of a factor of 3 and 2.5, 
respectively. 

Our results are consistent with dHvA measurements 
|l5| : while extended quasiparticles are seen deep into 
the vortex state in both NbSe2 and VaSi, the additional 
damping attributed to the superconducting gap below 
H C 2 increases more slowly in NbSe2- Corcoran et al. 
measured the electron-phonon constant \ e - p h — 0.3 for 
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